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Dielectric tensor and lattice dynamics of a-Pb02-type Ti02 have been investigated using the 
density functional perturbation theory, with a focus on responses of the vibrational frequencies 
to pressure. The calculated Raman spectra under different pressures are in good agreement with 
available experimental results and the symmetry assignments of the Raman peaks of Q;-Pb02-type 
Ti02 are given for the first time. In addition, we identified two anomalously IR-active soft phonon 
modes, Biu and Bsu, respectively, around 200 cm~^ which have not been observed in high pressure 
experiments. Comparison of the phonon dispersions at and 10 GPa reveals that softening of phonon 
modes also occurs for the zone-boundary modes. The Biu and Bsn modes play an important role 
in transformation from the Q;-Pb02-type phase to baddeleyite phase. The significant relaxations of 
the oxygen atoms from the Ti4 plane in the Ti2 02Ti2 complex of the baddeleyite phase are directly 
correlated to the oxygen displacements along the directions given by the eigenvectors of the soft Biu 
and B^u modes in the a-Pb02-type phase. 



I. INTRODUCTION 

Titanium dioxide (Ti02) has attracted much interest 
for a long time due to its potential applications in photo- 
catalysis. Besides the three well studied ambient phases 
(rutile, anatase, and brookite),^ ^ interest in Ti02 has 
also extended to high pressure phases to fully explore 
its novel structural and electronic propertiesj^.^ For in- 
stance, application of high pressure on Ti02 at room 
temperature induces a series of high-pressure phases with 
hardness approaching that of diamond^ The cotunnite- 
type Ti02 synthesized at pressure (P) above 60 GPa and 
temperature above 1000 K is believed to be the hardest 
oxide yet discovered^^ In addition, density functional the- 
ory calculation predicted that the high-pressure modifi- 
cation of Ti02 could lead to a polymorph which responds 
to visible light ^ extending the photocatalyst property of 
Ti02 from ultraviolet to visible light. 

The high pressure induced structural transition of 
Ti02 generally follows these sequences: rutile or anatase 
-> a-Pb02-type (space group: Pbcn) baddeleyite 
type (space group: P2i/c) orthorhombic structure 
(01, space group: Pbca) cotunnite type (Oil, space 
group: Prima) Critical pressures of transitions be- 
tween two successive phases are found to be strongly de- 
pendent on the size of the sample and the experimental 
temperature. ^^5^ Generally, at room temperature, the 
critical pressure of the transition from rutile (anatase) to 
a-Pb02-type phase is located in the region between 2.6 
and 10 GPa (2.5 and 7 GPa)^ii^-i^ Interestingly, it was 
found that all the high pressure phases recover to the a- 
Pb02-type phase rather than the initial rutile or anatase 
phase upon decompression to ambient conditions'^ 
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Since the a-Pb02-type phase plays a bridging role 
in connecting the ambient Ti02 phase and the higher 
pressure phases, a number of studies have been carried 
out on phase transformations involving the a-Pb02-type 
phase, such as rutile (anatase) ^ a-FhO^^^ and a- 
Pb02 ^ baddeleyite Among these, the mecha- 
nism for transition from the rutile phase to the a-Pb02- 
type phase has been well understood, and it can be de- 
scribed by the 1/2 [OlT] shear of every other hexagonal 
close-packed (hep) plane (Oil) of the rutile phase A 
soft phonon mode (Big) with a negative pressure deriva- 
tive was found in the vibrational frequencies^^ which is 
crucial to drive the initial structure to a CaCl2-type in- 
termediate structure during this transition j^^i^^ A similar 
mechanism was proposed by Sekiya for the anatase to a- 
Pb02 phase transitions^ In this case, a soft mode{Eg) 
observed around 200 cm~^ in the Raman spectra was sug- 
gested to be closely related to the phase transition. How- 
ever, for the subsequent a-Pb02 ^ baddeleyite phase 
transition, the understanding of which should form the 
basis of investigating the transitions of Ti02 to the higher 
pressure phases (01 and Oil), no phonon softening was 
observed, and the mechanism of this displacive phase 
transition has not been reported. 

In this study, we carry out density functional pertur- 
bation theory (DFPT) calculations to investigate the di- 
electric property and Born effective charge tensor of a- 
Pb02-type Ti02. In addition, we study its lattice dy- 
namical properties and their variations with pressure to 
identify soft phonon modes associated with pressure in- 
duced structural instability. The Raman intensity is cal- 
culated and presented for pressure from to 10 GPa and 
compared with the experimental spectrum. The sym- 
metry assignments of the Raman peaks for a-Pb02-type 
Ti02 are given for the first time. In contrast to the Ra- 
man active soft modes of rutile and anatase phases, two 
IR active soft phonon modes (^i^ and ^3^) with fre- 
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quencies around 200 cm~^ are identified for the a-Pb02- 
type phase. Based on the analysis of the two soft modes, 
we conclude that the TiaO triangular complex and the 
Ti202Ti2 complex are the relevant structural units for 
understanding the phase transformation of a-Pb02-type 
Ti02 ^ baddeleyite Ti02. 



II. METHODS 

All calculations were carried out using the CASTEP 
code^iiS^ with the Ceperley-Alder local density func- 
tional approximation.^^ Brillouin zone integration was 
performed using the Monkhorst-Pack A:-point sampling^! 
with 4x4x6 and 3x3x3 grids for a-Pb02- 
type and baddeleyite phases, respectively. Norm- 
conserving pseudopotentials^ were generated by the 
Opium package^ and employed throughout the study. 
For Ti atom, the pseudopotential was constructed to in- 
clude the 3s and 3p orbit als as semi-core states. The core 
radii Vc for the pseudo wavefunctions were tuned to be 
1.35 (35), 1.35 (3j9), 1.45 {3d), and 1.35 (4s) a.u. to en- 
sure a good transferability of the potentials. For O atom, 
the Tc values are 1.4 a.u. for all the 25, 2p and 3d va- 
lence states. The structures were relaxed until the forces 
exerted on the atoms are less than 0.01 eV/A, and the 
stress exerted on the cell is less than 0.02 GPa. 

The vibrational frequencies at the F point, the 
transverse-optical/longitudinal-optical (TO/LO) split- 
ting of the zone center optical modes, Raman spectra, 
and dielectric properties were calculated using the den- 
sity functional perturbation theory. Before calculations 
were carried out for the a-Pb02 and baddeleyite phases, 
phonon spectrum of r utile Ti02 was calculated and com- 
pared with available result^. It was found that for 
the chosen pseudopotentials, a plane- wave cutoff of 980 
eV was sufficient to obtain reliable frequencies of the F 
modes and this cutoff was used in subsequent calcula- 
tions. 

To obtain the phonon dispersion curves of a-Pb02- 
type Ti02, we first calculated the dynamical matrix 
based on a 4 x 4 x 4 grid of q points, with a (1/4,1/4,1/4) 
shift to include the F point, using the density functional 
perturbation theory. The Fourier interpolation scheme 
was used to generate the dynamic matrix at an arbitrary 
q point based on the aperiodic force constant matrix in 
real space produced by the cumulant sum method. The 
long-range dipole-dipole interaction was also taken into 
account by using the calculated anisotropic Born effective 
charge tensor and dielectric tensor. 



III. RESULTS 

A. Structure of a-Pb02-type Ti02 

Since both the vibrational properties and phase tran- 
sition of a compound are closely related to its atomic 
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FIG. 1: (Color online) Crystal structure of a-Pb02-type Ti02 
viewed from different directions to show (a) the TisO (filled 
triangles) and (b) Ti2 02Ti2 complexes (green lines), h de- 
notes the distance of the O atom from the plane defined by 
the Ti atoms. Red/small balls and blue/big balls denote oxy- 
gen and titanium atoms, respectively. 



structure, we present first calculated structural informa- 
tion of the a-Pb02-type phase. The relaxed lattice pa- 
rameters and atomic positions at zero pressure are given 
in Table I and are compared with available experimental 
values. The unit cell of a-Pb02-type Ti02 contains four 
formula units with four equivalent Ti atoms and eight 
equivalent O atoms. As shown in Fig. 1, each Ti atom 
is surrounded by 6 adjacent oxygen atoms with a coor- 
dination number (CN) of 6. The O atom is bonded with 
three Ti atoms forming an approximately "Y" -shaped 
TisO triangular complex, and two TiaO complexes fur- 
ther form a Ti202Ti2 complex (Fig. lb) through sharing 
a Ti2 edge, similar to the case of rutile Ti02. There are 
four Ti2 02Ti2 complexes in each unit cell. The centers 
of the four Ti2 62X12 complexes are located at a corner 
(0,0,0), the center (1/2,1/2,1/2), an edge (0,0,1/2), and a 
face center (1/2,1/2,0) of the cubic unit cell, respectively, 
which are the four inversion centers of the unit cell. In 
contrast to rutile phase where Ti and O atoms forming 
the Ti202Ti2 complex are all in the same plane, the TisO 
and Ti202Ti2 complexes in the a-Pb02 phase adopt a 
nonplanar structure. The two O atoms in the Ti202Ti2 
complex appear on opposite sides of the plane defined by 
the four Ti atoms (Ti4 plane). The formation of these 
nonplaner TisO complexes can be understood from the 
atomic relaxation driven by the Big soft phonon mode 
in the rutile phase under pressure ^2^*2^ These TisO and 
Ti202Ti2 complexes in the a-Pb02 phase play important 
roles in the analysis of Born effective charges and the de- 
scription of a-Pb02 ^ baddeleyite structural transfor- 
mation which will be discussed in the following sections. 

The pressure induced structural modifications, includ- 
ing lattice parameters and atomic positions are calcu- 
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FIG. 2: Variations of y/Vo and the average Ti-0 bond length 
of a-Pb02-type Ti02 with pressure. Here V is the volume at 
pressure P, and Vb is the volume at zero pressure. 



lated and results for P = 10 GPa are presented in Table 1. 
The distance (/i) of O atom from the Tis plane increases 
from 0.170 A at GPa to 0.195 A at 10 GPa. The varia- 
tions of y/Vb, where V is the volume at a given pressure 
and Vb is the volume at P = GPa, and the average 
Ti-0 bond length with pressure are shown in Fig. 2. A 
linear decrease are seen for both quantities. The bulk 
modulus {Bo) and its pressure derivative B are deter- 
mined through fitting of V/Vq to the Birch-Murnaghan 



equation 
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The fitted values of Bq and B are 248.6 GPa and 7.9, 
respectively. The bulk modulus obtained here is in good 
agreement with the experimental value of 258 ± 8 GPa^^ 



TABLE I: Calculated structural parameters, atomic fractional 
positions, and densities of a-Pb02-type Ti02 at and 10 
GPa. 



GPa 



10 GPa 



[a,b,c] [4.531, 5.465, 4.870] 

[4.531, 5.501, 4.906]^^ 
[4.61, 5.43, 4.87]2^ 
Ti (0, 0.176, 0.25) 

(0, 0.170, 0.25)2^ 
O (0.272, 0.384, 0.419) 

(0.272, 0.381, 0.414)^^ 
Density 4.401 
4.336^-^ 



[4.487, 5.374, 4.804] 

(0, 0.179, 0.25) 
(0.270, 0.389, 0.422) 
4.581 



B. Born effective charges and dielectric tensors of 
a-Pb02-type Ti02 

Born effective charge (BEG) tensor quantifies macro- 
scopic electric response of a crystal to the internal dis- 
placements of its atoms The calculated Born effective 
charges of the a-Pb02-type Ti02 are given in Table II. It 
is known that BEG can be strongly influenced by the local 
environment of each atom, and the form and the number 
of independent elements in the BEG tensors are deter- 
mined by the local atomic symmetry. Since a-Pb02-type 
possesses lower symmetry compared to rutile and anatase 
phases, the number of independent components in BEG 
tensors of Ti and O atoms in a-Pb02-type Ti02 is much 
larger than those of rutile and anatase phases ^^i^^ In ad- 
dition, in a-Pb02-type phase, since the Ti atoms occupy 
positions of higher local symmetry (C2) than O atoms 
(Ci), the number of independent elements in the BEG 
tensor for Ti atom is smaller than that of O atom. 

The principal directions (Pi,P2,P3) of the Born ef- 
fective charge tensor of the O atom are shown in 
Fig. 3, and the three principal values (^1,^2,^3) are 
(—5.35, —1.19, —3.01). The component in Pi direction 
is much smaller than the nominal ionic charge of the O 
ion (Z = —2). The Born effective charge tensor of the 
O atom is strongly anisotropic with the absolute values 
of Zi and Z3 much larger than that of Z2. This behav- 
ior arises from the bonding structure of O atom in the 
Ti30 complex. Each O atom in the Ti30 complex is 
bonded to three adjacent Ti atoms by forming sp^ -\ike 
bonds. Therefore, the polarizations along the Pi and P3 
directions are more sensitive to the displacement of O in 
the Ti3 plane, compared to that in the normal direction 
(P2). Since the Ti30 complex in the a-Pb02 phase has a 
similar structure as the Ti30 unit in the rutile case, the 
principal values of Born effective charge tensor for the O 
atom are comparable to the values (—5.06, —1.34, —3.96) 
of the O atom in rutile phase^^ 



TABLE II: Born effective charge, electronic (so 
(so) dielectric tensors of a-Pb02-type Ti02. 



and static 



Born effective charge tensors 



/ 6.89 0.14 \ 
Ti: 5.54 

\ 0.36 6.68 J 
/ -3.45 -1.28 -1.40 
O: -1.53 -2.77 
V -1.61 -0.07 -3.34 



Dielectric tensors 



£0 




The calculated electronic (£00) and static {eo) dielectric 
tensors are also given in Table II. The dielectric tensors 
are usually overestimated in local density approximation 
(LDA) calculations due to the underestimation of the 
band gap. However, application of the scissors correction 
has previously been found to improve the agreement with 
experimental data for both rutile and anatase casesj^ In 
this study, however, the dielectric tensors are calculated 
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FIG. 3: (Color online) Principal directions of the Born effec- 
tive charge tensor of the O atom in Q;-Pb02-type Ti02 (O 
atoms in red and Ti atoms in blue). The principal values 
along the three principal directions are —5.35, —1.19, —3.01, 
respectively. 

without the scissors correction. Compared with the pre- 
vious predictions*^ without scissor correction for rutile 
and anatase phase, we find that the electronic dielectric 
constant in the a-Pb02 phase is slightly smaller than 
those of rutile but larger than the values of the anatase 
phase, which may be attributed to a more similar struc- 
ture between a-Pb02 phase and rutile phase as shown in 
the Ti202Ti2 unit in both phases. 



C. Phonons of a-Pb02-type Ti02 

Since the primitive cell of a-Pb02-type Ti02 contains 
12 atoms, there are a total of 36 modes of vibration. 
The factor group analysis shows that the long wavelength 
optic phonon modes at the T point can be decomposed 
as 

All the phonon modes are non-degenerate. The Biu^ 
B2U1 Bsu modes are IR active, whereas the Ag^ Big^ B2g 
and Bsg modes are Raman active and the Au modes are 
silent. The Raman and IR modes are mutually exclusive 
due to the presence of inversion symmetry in the crystal. 

The calculated vibrational frequencies of the Raman 
modes, IR TO/LO modes and the silent Au modes, as 
well as their pressure dependence {doo/dP) are summa- 
rized in Table III. Here duj/dP is obtained through a 
linear fitting of uj — P data between and 10 GPa, which 
can be further used to generate the Griineisen parameter 
for each mode. An average value of 2.29 cm~^GPa~^ is 
obtained for duj/dP of the Raman modes. 

To enable a comparison with the experimental results, 
we calculated Raman spectrum of a-Pb02-type Ti02. 
The spectrum at GPa is shown in Fig. 4a. The cal- 
culated frequencies of Raman active modes in this work 
and the experimental values by Mammone et ali^ are 
indicated by the black and red lines at the bottom of the 
figure, respectively. By comparing the calculated and the 



TABLE III: Calculated phonon frequencies (oj) of Ti02 at 
the r point, and their pressure derivatives (duo/dP). For the 
IR-active modes, both TO and LO frequencies are listed. 



IR/Silent 


uo (TO/LO) 


du/dP 


Raman 




du/dP 










154.8 


0.67 


IR 






A. 


164.3 


0.13 


Biu 


197.4/229.9 


-2.61 


Big 


180.1 


1.71 


Bsu 


205.9/267.3 


-0.53 


Bsg 


215.3 


0.07 


Biu 


253.9/360.1 


3.62 


Big 


286.5 


-0.20 


Bsu 


280.5/286.3 


-0.06 


A. 


295.5 


-0.42 


B2u 


301.7/368.0 


0.54 


Big 


309.4 


1.18 


Bsu 


371.3/477.1 


5.97 


B2g 


344.6 


1.00 


B2U 


421.3/653.4 


2.74 


Bsg 


388.1 


2.67 


J^lu 


4Z ( .Z/ DZ ( .0 


Q 1 /I 
0.14 


^2g 


AC\A A 
4U4.D 


1.(0 


B3U 


483.1/838.7 


4.43 


A. 


431.9 


2.60 


Biu 


556.5/809.7 


3.78 


Big 


445.0 


4.53 


B2u 


701.0/788.5 


3.54 


B3g 


457.5 


2.60 








B3g 


154.8 


0.67 


Silent 






A. 


548.4 


5.06 


Au 


303.3 


-0.16 


Big 


576.3 


5.35 


Au 


527.2 


5.00 


B3g 


624.7 


4.69 


Au 


630.6 


4.00 


B2g 


813.0 


3.83 



experimental Raman spectra, we conclude that the fre- 
quencies at 150 and 172 cm~^ observed experiment all}*^^ 
can be assigned to the Bsg (154.8 cm~^) and Ag (164.3 
cm~^) modes, respectively. These assignments are fur- 
ther supported by the calculated high pressure Raman 
spectra shown in Fig. 4b. The gradual overlap of the po- 
sitions and the decrease in the intensities of the two peaks 
at high pressures are consistent with the trends observed 
in the experiment. Particularly, the disappearance of the 
peak (172 cm~^) observed in the experimental Raman 
spectra was used as an indicator of the occurrence of the 
phase transition from the a-Pb02-type Ti02 to badde- 
leyite phase. In addition, the splitting of the two peaks 
around 450 cm~^ {Ag and Big) and the significantly 
stiffening of three peaks above 500 cm~^ at high pres- 
sures were also clearly recorded in the experiments J^i^ 
It should be noted that, both the calculated and exper- 
imental Raman spectra show no soft modes with signif- 
icant negative derivatives of pressure, as the Big {Eg) 
mode in high pressure Raman spectrum of deformed ru- 
tile (anatase) phaseJ^*^ 

However, the two lowest IR active modes, Biu at 197.4 
cm~^ and B^u at 205.9 cm~^, respectively, show large 
negative pressure derivatives as seen in Table HI. Al- 
though several Raman and silent modes show similar 
negative-pressure derivatives, those of the two IR modes 
are much more stronger. These anomalously negative- 
pressure derivatives suggest a strong softening mecha- 
nism which can be correlated to atomic displacements 
in the structure under pressure and lattice instability. 
The frequencies of the two IR-active soft TO modes are 
shown in Fig. 5 as functions of pressure. A nonlinear 
dependence of the frequencies on pressure is clearly seen. 

The eigenvectors of the two soft modes Biu and Bsu 
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FIG. 4: (Color online) (a) Calculated Raman spectrum of a- 
Pb02-type Ti02 at GPa. The black and red lines at the 
bottom of the figure represent the calculated frequencies tab- 
ulated in Table III and the experimental results from Ref . [ij , 
respectively, (b) Raman spectra at various pressures from 
to 10 GPa. 
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FIG. 5: (Color online) Frequencies of the two IR-active Biu 
and Bsu soft modes as functions of pressure from to 10 GPa. 




(b) 



Bsu mode 




TI3O complex 1 TigO complex 2 
Ti202Ti2 complex 



FIG. 6: (Color online) Eigenvectors of the Biu (a) and Bsu 
(b) soft modes. The vibrational eigenvectors of the atoms in 
the TisO and Ti2 02Ti2 complexes (O atoms in red and Ti 
atoms in blue) are given on the right side of each figure. 



are shown in Fig. 6. The vibration patterns can be an- 
alyzed based on the Ti2 02Ti2 complex consisting of two 
TisO complexes through sharing a Ti2 edge as shown 
in the Fig. 1. The atomic vibrating patterns of the 
Ti30(Ti30-I and TisO-II) and Ti202Ti2 complexes are 
shown in the left part of the figure. The vertical move- 
ment of the oxygen atom out of the plane formed by the 
titanium atoms can be clearly observed. On the other 
hand, the oddness character (denoted by "ii" subscript) 
of Biu and Bsu under inversion requires that those atoms 
connecting with each other through the inversion symme- 
try operation have the same vibration components. Thus 
the two edge Ti atoms, connected with inversion opera- 
tion, vibrate in phase but slightly deviate with the vibrat- 
ing direction of the two end Ti atoms in the Ti202Ti2 
unit. Similarly, the two oxygen atoms locating on the 
two sides of the Ti4 plane also have the same vibrating 
component and they vibrate vertically to the Tia plane. 

To find out whether there exists any other soft phonon 
modes, we calculated the phonon dispersion curves and 
the vibration density of states (DOS) of the a-Pb02-type 
phase at GPa and 10 GPa, respectively. As shown in 
Fig. 7a, a significant shift of the spectrum towards the 
high frequency side is seen at frequencies above 300 cm~^ 
which is consistent with the larger pressure dependence 
of the high frequency phonon modes given in Table III. 
Fig. 7b shows details in the region between 150 and 250 
cm~^. These branches are related to the soft Biu and B^u 
modes. It is clear that in addition to the F point, fre- 
quencies of the modes along the T-Y^T-X ^T-Z and parts 
of Y-U^ Y-S^ X-S^ Z-T also decrease at high pressure. 
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FIG. 7: (Color online) (a) Calculated phonon dispersion 
curves of a-Pb02-type Ti02 at GPa (black solid lines) and 
10 GPa (red dashed lines). The phonon density of states 
(DOS) at the two pressure values are plotted in the left panel, 
(b) Details of the region from 150 to 250 cm~^ in (a). The 
vibrational branches associated with the IR active Biu and 
Bsu modes at 10 GPa are drawn in green lines to show mode 
softening. 




^ (II) 

FIG. 8: (Color online) Crystal structure of baddeleyite-type 
Ti02 (O atoms in red and Ti atoms in blue). The two types 
of Ti2 02Ti2 complexes in the structure are shown in brown 
color (type I) and green color (type II), respectively. The 
dashed lines depict the unit cell of the structure. The right 
panel shows the side views of Ti202Ti2-I and Ti202Ti2-II. 



D. Structure and Phonons of baddeleyite Ti02 

The atomic structure of the baddeleyite Ti02 (mon- 
oclinic) is illustrated in Fig. 8. The CN of Ti atom in 
this structure is 7, compared to 6 in a-Pb02-type Ti02. 
All Ti atoms are equivalent and occupy the 4e Wyckoff 
sites, but the 8 oxygen atoms can be divided into two in- 



equivalent groups, 01 and Oil, which occupy two sets of 
4e Wyckoff sites respectively. A structure similar to the 
Ti202Ti2 complex in a-Pb02 phase also exists in badde- 
leyite Ti02, as shown in Fig. 8. However, in baddeleyite 
Ti02, there are two types of Ti202Ti2 complexes, termed 
as Ti2 02Ti2-I and Ti2 02Ti2-II, respectively, depending 
on the type of oxygen (01 or Oil) involved. The distances 
from the oxygen atoms to the Ti4 plane in the two types 
of Ti202Ti2 complexes are different. Using the experi- 
mental lattice parameters of a = 4.640 A, 6 = 4.760 A, 
and c = 4.810 A^ the 01 and Oil atoms in the relaxed 
structures are 0.358 A and 0.940 A from the Ti4 plane, 
both being much larger than the corresponding value of 
0.12 A in the a-Pb02-type Ti02 at zero pressure. 

Finally, the calculated phonon frequencies at the F 
point of baddeleyite Ti02 are listed in Table IV. Fif- 
teen IR active {SAu-\-7Bu) and eighteen Raman active 
{9Ag-\-9Bg) modes are expected for baddeleyite phase 
from the factor group analysis. 

TABLE IV: Calculated phonon frequencies of IR active modes 
{Au and Bu) and Raman active modes {Ag and Bg) for bad- 
deleyite Ti02 at the F point. For the IR- active modes, both 
TO and LO frequencies are listed. 

IR/Raman Frequency (cm~^) 

Au (TO/LO) 227.9/242.8, 246.1/256.7, 331.9/332.4, 

384.1/404.5, 439.7/540.5, 451.6/452.3, 

582.9/860.7, 733.3/736.3 
Bu (TO/LO) 110.5/196.3, 295.8/341.1, 359.2/457.9, 

495.7/520.0, 529.9/560.8, 561.5/702.5, 

784.7/960.0 

Ag 136.8, 217.0, 249.7, 348.7, 399.4, 427.0, 

511.7, 647.3, 678.5 
Bg 189.1, 296.6, 352.8, 423.2, 430.4, 493.5, 

621.7, 707.0, 814.9 



IV. DISCUSSION 

The transition from the a-Pb02-type phase to the bad- 
deleyite phase mediated through the soft phonon modes 
can be derived directly from the comparison of the re- 
lated softening modes and the structural features of bad- 
deleyite phase. The significant deviation of the oxygen 
atoms from the Ti4 plane in the Ti202Ti2 complex of 
the baddeleyite phase arises from relaxation of oxygen 
atoms along the directions described by the eigenvectors 
of the Biu and Bsu soft modes in a-Pb02 Ti02. This soft 
phonon mode based mechanism was further confirmed by 
results of our molecular dynamics (MD) simulation of the 
a-Pb02-type Ti02 which was carried out at experimen- 
tal conditions of 300 K and 37 GPa.^^ The MD simu- 
lation was implemented by the "General Utility Lattice 
Program" (GULP) package^^ with a NPT ensemble using 
Dreiding empirical force potential^ which has previously 
been adopted in simulations of molecular adsorption on 
Ti02 surface, zeolite confined Ti(0H)4 nanoparticles^ 
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FIG. 9: (Color online) Comparison of the structure of a- 
Pb02-type Ti02(a) with that of the deformed Q;-Pb02-type 
Ti02(b) obtained from molecular dynamics simulations. The 
right panel shows the side views of the structures of Ti202Ti2. 
Red (small) balls and blue(big) balls represent oxygen and ti- 
tanium atoms, respectively. 

and other surfactant interactions at mineral surfaces 4^ 
The Nose-Hoover thermostat was applied to control the 
temperature. It should be noted that the initial veloc- 
ities, which were produced randomly with probabilities 
based on Gaussian distribution, are not related to the 
mode direction of Biu and B^^u soft modes. The ini- 
tial structure and the snapshot of a representative in- 
termediate structure of a-Pb02-type Ti02 are shown in 
Fig. 9^ Strong distortions in the lattice and out-of-plane 
displacement of the oxygen atoms can be observed. The 
single type of Ti2 62X12 complex in the a-Pb02 phase 
evolves into two types (I and II) which are similar to 
those in the baddeleyite structure given in Fig. 8. The 
relaxation of the Ti202Ti2 complex which shows charac- 
teristic features of the Biu or B^^u soft phonon induces 
significant shear deformation in the y-z plane. It should 
be pointed out that a complete structural transformation 
from the a-Pb02 phase to the baddeleyite phase was not 
observed in this simulation. This is due to constraints of 
the supercell (2x2x2) and possible high potential barrier 
along the pathway. Nevertheless, result of the molecular 
dynamics simulation does indicate that the local atomic 
displacements are closely related to the Biu and B^^u soft 
modes at the F point. 

In this study, the structural difference between the 
a-Pb02-type and the baddeleyite phases is analyzed 
based on the basic building units TiaO and Ti202Ti2 
complexes, in contrast to the commonly used TiOe 
octahedra4^i^ One clear advantage of using TiaO and 
Ti202Ti2 as the basic units is in the description of the 



Biu and B^u soft modes of the a-Pb02-type phase, and 
the correlation between the phase transition and the out- 
of-plane displacements of the oxygen atoms in these soft 
modes. As a matter of fact, such soft modes not only 
dictate phase transition initiated from the a-Pb02-type 
phase, similar features have also been observed in the 
high pressure phase transitions originating from rutile 
and anatase phases of Ti02. The Big soft mode in ru- 
tile phase and the Eg soft mode in anatase phase can be 
described similarly in terms of TiaO and Ti202Ti2 com- 
plexes, where out-of-plane displacements of the oxygen 
atoms relative to the Tia plane are observed. ^^^^^ For in- 
stance, in the case of rutile to a-Pb02 phase transition, 
the displacements of the O atoms associated with the 
Big mode lead to the CaCl2 structure with corrugated 
Ti202Ti2 unit chains along the [110] direction, and the 
final a-Pb02-type structure can be achieved by a sub- 
sequently shearing in the (Oil) plane of the CaCl2-type 
structure In the anatase case, the high pressure 
transition to the a-Pb02 phase has also been found to 
be closely related to the soft Eg mode, and the gliding of 
anionic layers results in the transformation.— Thus, the 
instability of the TiaO units may be regarded as a com- 
mon feature of high pressure transformations of Ti02. 



V. CONCLUSIONS 

Softening of phonon modes, which is usually an impor- 
tant indicator of structural instability during displacive 
phase transition, has previously been found in both ru- 
tile and anatase Ti02 under high pressure. However, 
for phase transitions originating from the a-Pb02-type 
Ti02, no soft modes have yet been reported during high 
pressure transition. Using the state-of-art DFPT calcula- 
tions, we show that there exist two IR- active soft phonon 
modes {Biu and B^u) around 200 cm~^ which are absent 
in the high pressure Raman spectrum. These two modes 
play an important role in understanding the atomic dis- 
placements in a-Pb02-type Ti02 and transformation to 
baddeleyite polymorph. In addition to zone-center fre- 
quencies, the changes in phonon dispersions with pres- 
sure indicate that significant softening also occurs in the 
zone-boundary modes, which suggests that gliding of lay- 
ers in a-Pb02-type Ti02 is likely to occur during the 
a-Pb02-type Ti02 phase transition. Furthermore, the 
symmetry assignments of the Raman peaks for a-Pb02- 
type Ti02 are given for the first time. Since there is usu- 
ally a coexistence of different phases during the process 
of high pressure phase transitions, the calculated Raman 
spectrum and its change with pressure are valuable for 
experimentalist to identify the onset of phase transition 
related to the a-Pb02-type phase. 
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